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The number of pseudoscalar mesons in the mass range from 1400 to 1500 MeV/c^ has been a 
^ . subject of considerable interest for many years, with several experiments having presented evidence 

Cn ' for two closely spaced states. A new measurement of the reaction 7r~p — > K'^ K'n'^n has been made 

J^^ \ at a beam energy of 18 GeV. A partial wave analysis of the K^ K^-k*^ system shows evidence for 

three pseudoscalar resonances, ?7(1295), r;(1416), and r;(1485), as well as two axial vectors, /i(1285), 
and /i(1420). Their observed masses, widths and decay properties are reported. No signal was 
observed for C(1480), an I J^ = 1^1 state previously reported in 07r" decay. 
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A long history of confusion surrounds the pseudoscalar 
meson spectrum |n^,p| . Results from previous experiments 
|p|-|l8[ carried out over the past four decades, including 
peripheral production with pion beams, pp annihilation, 
and J/ip radiative decay, indicate that there may be three 
isosinglet pseudoscalar states in the narrow mass range 
from 1250 to 1500 MeV/c^. However only two states are 
expected in the standard J^ = 0~ nonet. 

The lowest pseudoscalar, 77(1295), was first observed 
in a PWA analysis of the 'qmr system [091 and later 
confirmed by other experiments p0|-|2^. However the 
77(1295) was not observed in some reactions where the 
production of pseudoscalars is allowed, such as 77 col- 
lisions IgJ], central production [g5 26 1 and J/ip decays 
||10|-[l2|. Therefore further confirmation of this state is 
desirable. 

The central issue in the higher mass region is some- 
times referred to as the E/l puzzle. There was a long- 
standing question whether the iota, seen in radiative J/tp 
decay ||l|,0, and the E meson (now called /i(1420)) 
were the same state or two distinct resonances. Pe- 
ripheral production experiments at Brookhaven National 
Laboratory showed a complex spectroscopy in this mass 
region |ij2^,|§. 

One of these, experiment E771, used pion, kaon and 
antiproton beams to study the production of different 
states [|J2|,||]. The kaon data showed that the /i(1285) 
has very little strange-quark content pS]. Evidence was 
also obtained for a J^'~'' = 1"' resonance at about 1380 
MeV/c^ [|8|, below the K*l< threshold. Those resuhs 
were consistent with data obtained by the LASS collabo- 
ration ||2^. They also observed J^^ = 1++ states at 1420 
and 1500 MeV/c^ 0, and two pseudoscalar states [|^J28| 
in the 1400 MeV/c mass region, one decaying mainly to 
K*K and the other decaying primarily to ao(980)7r. Ev- 
idence was also presented for an additional pseudoscalar 
at 1515 MeV/c^ appearing only in ao(980)7r 0. 

In their latest publication ||ll[ the Mark HI group con- 
clude that their data require two J^ = 0^ resonances: 
one at 1416 MeV/c^ decaying to ao(980)7r and the other 
at 1490 MeV/c^ decaying to K*K. The DM2 collabora- 
tion [|2| present quite a different picture: a 0~ resonance 
at 1421 MeV/c^ decaying to K*K and a 0~ resonance at 



1459 decaying to ao(980)7r. The Obelix group recently 
reported results for pp -^ {K'^KsTr^)7T^n'' at rest [0-^. 
They observe two 0~ resonances in the (KKtt) final state, 
one at 1405 MeV/c^ decaying directly to (KKtt) and 
a broader one at 1500 MeV/c^ decaying primarily to 
K*K. The PDG concludes that most of the evidence 
now supports two 0^ states in the range 1.4-1.5 GeV/c^, 
but their decay properties are poorly determined. The 
present experiment was designed to clarify this situa- 
tion by obtaining a large sample of data for the reaction 
TT^p -^ K^K^Ti^n. The task is a challenging one be- 
cause of the large number of overlapping resonances in 
this mass region, and the presence of the K*K threshold 



at about 1.4 GeV/c^. 

One controversial state is the C(1480), which was 
found by the Lepton-F group in (/)(1020)7r'' decay [ p0[ . 
The mass and width of this resonance were measured to 
be M = 1480±40 MeV/c^, and F = 130±60 MeV/c^, re- 
spectively, with I'-^J^'-' = 1+1^^. Because it has a large 
production cross section in Tr~p -^ {K~^K~Tr°)n and its 
decay into cutt^ was not observed, it has been suggested 
that this state is not a conventional meson [pil. The 
PDG now lists C(1480) under /9(1450), a state previously 
associated with non-strange decay modes [El. 

The present experiment (E852) was carried out dur- 
ing the 1997 running period of the Alternating Gradient 
Synchrotron at Brookhaven National Laboratory. A di- 
agram of the experimental apparatus is shown in Fig. 1. 
A Cherenkov-tagged 7r~ beam of momentum 18.3 GcV/c 
was incident on a 60 cm liquid hydrogen target within the 
Multi-Particle Spectrometer. The spectrometer central 
field was set to 0.75 Tesla. The target was surrounded by 
a four-layer cylindrical wire chamber (TCYL) which was 
used to reject events that included a large-angle charged 
particle |g^. A segmented Csl detector was used to re- 
ject events with large-angle photons | p3[ . 

The downstream part of the magnet was equipped 
with six multi-plane drift-chamber modules (DCl-6) |3J] 
for charged-particle tracking. Three proportional wire 
chambers (TPXl-3) allowed the acquisition system to 
trigger on events having two charged tracks. A 3045- 
element lead-glass calorimeter (LGD) pa] was used to 
detect forward photons, and a 96-segment threshold 
Cherenkov counter was used to distinguish charged kaons 
from pions. The Cherenkov radiator was Freon 114 at at- 
mospheric pressure, giving a diffractive index of 1.00153. 
The corresponding pion (kaon) threshold is at 2.52 (8.93) 
GeV/c. 

The trigger for the present data required two forward- 
going charged tracks, energy deposition in the LGD, and 
no charged recoil track. A total of 18.5 million events 
of this type were recorded. After calibration and track 
fitting the events were filtered by requiring: 

1. a fully reconstructed beam track; 

2. two forward tracks, one positive and one negative, 
both identified as kaons by the Cherenkov detector; 

3. two energetic clusters in the active region of the 
LGD, both identified as neutral particles; 

4. a reaction vertex within the target volume; 

5. energy deposition < 40 MeV in the Csl detector. 

In the final step of the data selection a kinematic fit 
was made to the tt~p — > K^K~n^n reaction hypothe- 
sis using the full covariance matrix from track and ver- 
tex reconstruction [Bm . Events with confidence level less 
than six percent were rejected. This was done to remove 



events with low-energy photons that were not detected 
by the apparatus. Approximately 34,000 events survived 
the above selection criteria. Fig. 2 shows the invari- 
ant mass distributions for these exclusive events. The 
most striking features of the data are two peaks in the 
K~^K~TT° mass spectrum at about 1.3 and 1.45 GeV/c^. 
The two-particle mass spectra show clear signatures for 
0(1020) and K*{892) excitation. 

A partial wave analysis (PWA) of the present data was 
made in the isobar model i3^,pql . The mass of the three- 
meson final state was binned in 20 MeV/c^ intervals and 
independent fits were performed on the data in each bin. 
The final state was represented as a sequence of interfer- 
ing two-body intermediate states. An initial decay of a 
parent meson into a KK or Ktt'^ intermediate resonance 
(isobar) and an unpaired pion or kaon, followed by the 
subsequent decay of the isobar, populates the K^K~Tr° 
final state. Each partial wave is characterized by the 
quantum numbers J^'~' [isobar]LM'^ , where J^ are the 
spin, parity and C-parity of the partial wave, M is the 
absolute value of the spin projection on the beam axis, e 
is the reflectivity (corresponding to the naturality of the 
exchanged particle), and L is the orbital angular momen- 
tum between the isobar and the unpaired particle. Both 
nucleon spin-flip and non-flip amplitudes were included. 
The experimental acceptance was determined by 
means of a Monte Carlo simulation, which was then in- 
corporated into the PWA normalization for each partial 
wave. In the description of the isobars a coupled-channel 
Flatte parametrization p9] was used for the ao(980). The 
parameters of this function were set to the values used in 
previous partial wave analyses 0,0] . A scattering length 
function was used to describe {Kt:)s, the (Kn) S-wave in- 
teraction 0. The other isobars, K*{892), and 0(1020), 
were modeled by relativistic Breit-Wigner shapes with 
masses extracted from the PDG [Q. Experimental res- 
olution made a significant contribution to the measured 
width of the 0(1020) so in that case the observed width 
(9.5 MeV/c^) was used in the PWA analysis. 

In the partial wave analysis numerous fits with differ- 
ent wave sets were performed to determine the minimum 
set of waves that gave a good description of the angu- 
lar distributions and the two-particle mass distributions. 
All waves were included as C-parity eigenstates ^^ (eg. 
X*+K^ ± K*^K+). Partial waves with J < 4 were in- 
cluded for K* , gq and isobars as well as {Ktt)s- The 
extended maximum likelihood method was used to de- 
termine the goodness-of-fit. 

The largest contributions came from pseudoscalar and 
axial vector waves, with natural-parity exchange domi- 
nating the results. Table I lists all partial waves used 
in the final fit. Waves with J^^ = 2+"*" and 1~~ were 
employed to account for the presence of strong tails from 
the /3(1700) and other high- lying states. Both of these 
were needed in the fits. The 1"' K*K S-wave was in- 
cluded to allow for possible excitation of the /ii(1380). A 



non-interfering isotropic background wave was included 
at each mass bin to simulate the cumulative effect of nu- 
merous small waves that were omitted from the fit. 

The final fit was performed on 20,000 events having 
momentum transfer \t\ > 0.1 GeV^. This cut reduced 
the contributions from 07r decay (see Fig. 2), and also 
the strength of the 1 waves. 

No {K'k)sK waves were included in the final fit. Those 
waves were found to be ambiguous with the oott S-waves 
and with the background wave. Fits which included both 
{Kt:)s and ao waves exhibited large interference that var- 
ied wildly in strength at adjacent M{KKn) mass bins. 
Therefore one cannot exclude the possibility of some un- 
detected {KTr)sK strength in the present analysis. 

Different wave sets were used to describe the data 
above and below the K*K threshold. Fits below 1.375 
GeV/c^ did not include K*K waves. The main results 
of the mass- independent PWA fit are shown in Figs. 3 
and 4. The low-mass spectrum shows large contributions 
from J^'^ = l"*""*" and '" waves. 

Mass-dependent fits to the intensities and phase differ- 
ences were made by least-squares minimization for three 
separate wave pairs. The phase information is partic- 
ularly important for determining the resonant content 
of the mass-independent results. These mass-dependent 
fits used linear combinations of relativistic Breit-Wigner 
poles with mass-dependent widths and Blatt-Weisskopf 
barrier factors |El|. The details of this procedure are 
given below. 

The first fit was to the oott" waves below 1.36 GeV/c^. 
Those results, depicted in Fig. 3 (a-c), show clear exci- 
tation of r/(1295) and /i(1285). The resonance param- 
eters from this fit are tabulated in Table II. Contribu- 
tions from experimental resolution (about 10 MeV/c^) 
were removed from the listed widths. The errors given in 
the table include statistical contributions from the mass- 
dependent fit (first entry), and systematic errors from 
the PWA fits (second entry). The latter contribution 
was estimated by comparing the results of several mass- 
dependent fits which used different PWA results. In those 
fits the mass binning and wave sets were varied. 

With the exception of the /i(1285) width, the mea- 
sured parameters of the ry(1295) and /i(1285) are in good 
agreement with the average values reported by the PDG 
§. For the /i(1285) we observe a width F = 45 ± 9 ± 7 
MeV/c^, which is consistent with the latest measure- 
ment from GAMS |E3| but larger than the PDG value 
of 24.0 ± 1.2 MeV/c\ Fitting only the intensity func- 
tion in this case (data in Fig. 3b) yields a much smaller 
width, r = 23 ± 5 MeV/c^ Thus it is the interference 
with 77(1295) that is important in our determination of 
the width. 

Above 1.38 GeV/c^ several overlapping states show in- 
teresting structure, as indicated in Fig. 4. Pseudoscalar 
waves show up prominently in aoTr" and K* K, while the 
most important axial vector wave is in K* K decay. The 



pseudoscalar mass spectrum for uott^ decay is markedly 
different from that for K*K. The former shows a narrow 
peak in the intensity distribution while the latter shows 
strength over a much wider mass range. 

The results of the mass-independent PWA fits for this 
mass region were interpreted in two steps. In the first 
step the intensities and relative phase of the most promi- 
nent pseudoscalar ao7r° and axial vector K*K waves were 
fitted with two Breit-Wigner poles. The results of that 
fit are given in Fig. 4 (a-c) and in Table II. They are 
consistent with previous observations of /i(1420) and a 
low-mass component of 77(1440) H, seen here at a mass 
of 1416±4±2 MeV/c^. In the present work this state will 
be labeled r;(1416). The mass and width of the /i(1420) 
agree with the average values reported by the PDG 0. 

In the second step of the mass-dependent analysis 
the intensities and relative phase of the prominent pseu- 
doscalar and axial vector K*K waves were fitted to three 
resonance poles. The parameters of the /i(1420) and 
77(1416) were held constant at the values determined 
above (see Table II) , while the mass and width of a third 
resonance were varied to minimize x^ ■ The results of this 
fit are given in Fig. 4 (d-f) and in Table II. They show a 
mass and width for the high-mass component of 7/(1440) 
equal to 1485±8 ± 5 MeV/c^ and 98±18 ± 3 MeV/c^, re- 
spectively. This state will be labeled 77(1485). For this fit 
X^ /d.o.f was 1.07 (for 16 degrees of freedom). Fits were 
also made with the parameters of the lower states uncon- 
strained. Those results were in good agreement with the 
values given in Table II, but with larger error bars. Pre- 
vious PWA analyses have determined that the K^K tt" 
spectrum is dominated by positive G-parity excitation 
iQ,^, so one is justified in assigning isospin / = to the 
positive C-parity states discussed above. 

The present results agree with those from earlier 7r~p 
[g[, pp 1^, and J/ij] radiative decay ||ll[ experiments; two 
pseudoscalar resonances were observed, with the lower 
state appearing clearly in ai^ir^ decay and the upper one 
strongly populating only the K*K spectrum. In the 
present measurements 7/(1416) was observed in both aoTr" 
and K*K decay. In order to test the significance of this 
result a new fit of the intensities and relative phase of the 
waves depicted in Fig. 4 was made. In that fit a single 
resonance was assumed in the pseudoscalar K*K wave. 
This resulted in the width of the resonance increasing 
to about 300 MeV/c^, and y^ /d.o.f increasing from 1.07 
to 1.38 (for 17 degrees of freedom), which demonstrates 
that the two-pole fit offers an improved description of the 
data. 

No strong signals were observed for excitation of the 
/i(1510). As one can see in Fig. 4, the J'"^ ^ 1++ K*K 
intensity shows a small bump at 1.5 GeV/c^ but no clear 
phase motion with respect to the 7/(1485). 

The decay properties of the 7/(1416) are of particular 
interest for determining the structure of the state. From 



our measurements we obtain the ratio of branching frac- 
tions, 



Bi?[77(1416) -^ K*K,K* -> Kn^] 
BR[riiU16) -> aoTrO, ao -> K+R- 



= 0.4 ±0.1 



When all charge states are included this leads to the ratio 
of total branching fractions, 



Bi?[77(1416) ^ K*K + c.c.] 
Si?[7/(1416) -^ ao7r] 



= 0.084 ±0.024 



The quoted error is statistical in nature. This ratio 
is reduced by about a factor of two when pseudoscalar 
{K-k)sK waves are included in the fits. The above val- 
ues are consistent with early results obtained by Obelix 
when similar wave sets were used [^ , but is much smaller 
than those obtained in later Obelix experiments [0,|| . 

The 7/(1416) and 77(1485) states were both observed 
to decay to K*K. The ratio of production rate times 
branching fraction for these two states is, 



i?[77(1416) -> K*K] 
i?[77(1485) -^ K*K] 



=4 =0.16 ±0.04 



Previous tt p experiments were unable to detect the 
?7(1416) contribution to the K*'K spectrum ||,|,|7)- 

In principle one should be able to use the measured 
widths and decay rates of 77(1416) and 77(1485) to deter- 
mine their valence structure. Unfortunately the present 
state of hadron theory does not allow this. In the "^Pq 
model one expects a pseudoscalar ss meson to have a 
narrow width (< 100 MeV/c^) and a dominant K*K de- 
cay branch 1421. However both 7/(1416) and 7/(1485) have 
widths that fall within acceptable limits, and both decay 
to K*K. The fact that 7/(1485) decays primarily by kaon 
emission suggests that it is the better ss candidate, but 
further theoretical work is needed before an identification 
can be made. 

Another important spectroscopic issue is the existence 
of C(1480). That state was first identified in a pro- 
duction experiment similar to the present one pG]. It 
was observed in 0(lO2O)7r° decay with a mass and width 
M = 1480 ± 40 MeV/c2, and F = 130 ± 60 MeV/c^ 

Subsequently it was not seen in several other experi- 
ments, including pp central production P5[ and pp an- 
nihilation at rest ^^. A 0(1020) side-band study and 
partial wave analysis were used to search for such a state 
in the present data. Both studies were carried out using 
the full data set without a t cut. 

In the K'^K~ effective mass spectrum a clear peak is 
seen at the mass of the </> meson (Fig. 2b). The intensity 
peaks at 1020.3 ± 0.4 MeV/c^ and has a width of 9.5 ± 
1.0 MeV/c^. The width is mainly from the instrumental 
resolution of the detector. For the side-band study the 
two- kaon mass spectrum was divided into three bins, each 
10 MeV/c^ wide, and K^K~tt^ mass distributions were 



produced for each bin. Those spectra are shown in Fig. 
5. Fig. 5b shows the K~^K~tt'^ effective mass spectrum 
for those events in the (f> signal region. Background was 
removed from this spectrum by subtracting the average 
number of events in the adjacent two-kaon mass intervals. 
The resulting background-subtracted mass spectrum is 
shown in Fig. 5d. No peaks can be discerned in the 
spectrum. 

A more definitive test was made by performing addi- 
tional PWA fits with (/)7r° S and P waves added to the 
wave set given in Table I. No peaks were observed in the 
resulting mass spectra. Since the C(1480) was believed 
to be a vector meson, these fits included all allowed ikf^ 
values for 1 (fin'^. Fig. 6 shows the contribution from 
the wave with the largest intensity (Af^ = 0^). 

The present results indicate no strong excitation of 
any (J^tt" waves and negligible resonance strength for 
1 07r° waves in particular. The previous identification 
of C(1480) cannot be confirmed. 

In summary, a partial wave analysis of the mesons pro- 
duced in the reaction n^p -^ K^K^n'^n was performed. 
The intensity and phase of the resulting waves show 
clear excitation of several previously identified states: 
/i(1285), ?7(1295), ?7(1416), /i(1420), and r/(1485). The 
existence of three low-mass pseudoscalars is confirmed. 
The ?7(1416) and ?7(1485) are distinguished by their de- 
cay properties; ?7(1416) decays primarily to aoTr" but has 
a small K*K branch as well. The 77(1485) was observed 
only in K*K decay. C(1480), a state previously identi- 
fied in KKtt decay, was not observed. 
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FIG. 1. Schematic diagram of the experimental apparatus. 
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FIG. 2. Mass spectra for (a) K+K-tv°, (b) Tf+i^", (c) 
K^tt'^, and (d) /(""tt". The open areas are all events and the 
shaded areas are events with \t\ > 0.1 GeV^. 



FIG. 3. Results of mass-independent fits (points) and 
mass-dependent fits (lines) for (a) '' [ao(980)] SO"*" inten- 
sity, (b) 1^"*" [ao(980)] PO'*' intensity, (c) phase difference be- 
tween (a) and (b). 
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FIG. 4. Results of mass-independent fits (points) and 
mass-dependent fits (lines) for (a) '" [ao(980)] SO""" inten- 
sity, (b) 1++ [A'* (892)] S0+ intensity, and (c) phase difference 
between (a) and (b), (d) 0"+ [A'* (892)] P0+ intensity, (e) 1++ 
[K* (892)] SO'^ intensity, and (f ) phase difference between (d) 
and (e). 
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TABLE II. Resonance parameters and decay modes of the 
observed states. Statistical errors are listed first, followed by 
systematic errors. 



Resonance 


M (MeV/c^) 


r (MeV/c^) 


Decay Modes 


/i(1285) 


1288 ± 4 ± 5 


45 ± 9 ± 7 


ao7r" 


r,(1295) 


1302 ± 9 ± 8 


57 ±23 ±21 


aoTT 


r,(1416) 


1416 ±4 ±2 


42 ± 10 ± 9 


aoTT°,K''K 


/i(1420) 


1428 ± 4 ± 2 


38 ± 9 ± 6 


K*K 


r;(1485) 


1485 ± 8 ± 5 


98 ±18 ±3 
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FIG. 5. 0(1020) side-band study. The K+R-'k'-^ effec- 
tive mass for events with (a) 1.0 < M{K+K') < 1.01 
GeV/c^ (b) 1.015 < M{K+K-) < 1.025 GeV/c^ (c) 
1.03 < M{K+K-) < 1.04 GeV/c^ and (d) the back- 
ground-subtracted signal. 
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FIG. 6. Intensity of the 1 — [<;/)(1020)]P0" wave. 



TABLE I. Partial waves used in the amplitude analysis. 
Note that the K* K partial waves were used only for masses 
greater than 1.375 GeV/c^ 
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